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ABSTRACT 



The mass density field in the local universe, recovered by the potent method from 
peculiar velocities of ^3000 galaxies, is compared with the density held of optically 
selected galaxies. Both density fields are smoothed with a Gaussian filter of radius 12ft, -1 
Mpc. Under the assumptions of gravitational instability and a linear biasing parameter 
bo between optical galaxies and mass, we obtain (3 = ^ 06 /b o = 0.74±0.13. This result 
is obtained from a regression of potent mass density on optical density after correcting 
the mass density field for systematic biases in the velocity data and potent method. 
The error quoted is just the la formal error estimated from the observed scatter in the 
density-density scatterplot; it does not include the uncertainty due to cosmic scatter 
in the mean density or in the biasing relation. We do not attempt a formal analysis of 
the goodness of fit, but the scatter about the fit is consistent with our estimates of the 
uncertainties. 
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1 INTRODUCTION 

With the advent of large-scale redshift and peculiar veloc- 
ity surveys, it has become possible to investigate the rela- 
tionship between the distribution of mass and galaxies on 
supercluster scales. Assuming that structures grow through 
gravitational instability (GI), then, on large scales where lin- 
ear theory is applicable, the peculiar velocity of a galaxy is 
proportional to its peculiar gravitational acceleration. The 
constant of proportionality depends on the density parame- 
ter, fi, like /(fi) ~ n 6 (e.g. Peebles 1980). The acceleration 
scales with the amplitude of the mass density fluctuation 
field, so, if we assume that the density contrast of optically 
selected galaxies, So, is related to the density contrast of 
mass, 8, by linear biasing So = boS, then the appropriate 
scaling factor between peculiar velocity and galaxy density 
becomes j3o = f(Q,)/bo- The determination of f3o is the goal 
of this paper. 

There are two general ways to obtain f3o from pecu- 
liar velocity data and the density field as derived from a 
galaxy redshift survey. Each method has advantages and 
disadvantages. One approach is to predict radial peculiar 
velocities from the density field via gravity and adjust f3o 
for the best fit to the peculiar velocity data (e.g. Huchra 
1988 and references therein; Strauss 1989; Kaiser et al. 1991; 



Shaya, Tully & Pierce 1992; Hudson 1994b; Willick et al., 
in preparation). A potential error in such velocity-velocity 
comparisons arises from the non-locality of the predicted ve- 
locities, which are affected by the density field in unsurveyed 
regions such as the Galactic Zone of Avoidance and distant 
regions. Another potential source of error is the fact that the 
computation of the predicted velocities involves smoothing, 
while they are compared to raw, unsmoothed peculiar ve- 
locity data. 

An alternative approach is to compare density fields 
smoothed in a similar way. In linear theory, the mass density 
fluctuation can be determined from the divergence of the 
peculiar velocity field, v, 



(i) 



and a generalized approximation is valid in the quasi-linear 
regime (see Section ^| below) . The comparison of densities in 
the linear regime yields j3o, and in the quasi-linear regime it 
allows a determination of 6o for a given value of Q. The main 
advantages of this approach are the locality of the compari- 
son (on scales larger than the smoothing scale) , and the fact 
that similar smoothing is applied to both fields at a scale 
on which linear or quasi-linear theory is applicable. This ap- 
proach was taken by Dekel et al. (1993, hereafter DBY93), 
who compared the mass density derived by potent (Dekel, 
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Bertschinger & Faber 1990, hereafter DBF; Bertschinger et 
al. 1990) from the peculiar velocity sample of Burstein (1990, 
Mark II) to the density field of galaxies drawn from the IRAS 
1.9- Jy redshift survey (Strauss et al. 1990, 1992). Dekel et 
al. used an elaborate likelihood analysis to find that the 
data were consistent with the assumptions of GI and linear 
biasing, and estimated fii = 1.3 ± 0.3. 

The goal of this paper is to estimate fio, under the 
assumption of GI and linear biasing, by a straightforward 
comparison of the density field of optically selected galaxies 
and the mass density field as determined by an improved 
potent procedure applied to an extended and newly cali- 
brated data set. 

The outline of this paper is as follows. The mass density 
field is first determined by applying potent to a data set 
of peculiar velocities containing ~3000 galaxies. The veloc- 
ity data, some details of the potent method, and in par- 
ticular the biases introduced in using sparse, non-uniform 
and noisy radial peculiar velocities to obtain a mass den- 
sity field are summarized in Section ^. The optical density 
field is then described in Section ^. In Section [|, we estimate 
Po by means of a linear regression of the potent density 
fluctuation, Sp, on the optical density fluctuation, So, elim- 
inating biases in the potent procedure, and evaluate the 
robustness of the result. We also demonstrate the difficulty 
of separating SI and feo via non-linear effects. Finally, in Sec- 
tion f| we summarize our results and discuss them in light 
of other determinations. 



2 POTENT DENSITY FIELD 

The sample of galaxies with measured redshifts and forward 
Tully-Fisher or D n — a (hereafter collectively TF) inferred 
distances used in this paper is a preliminary version of the 
Mark III catalogue. The calibration and compilation of the 
several data sets comprising the Mark III catalogue are de- 
scribed in detail by Willick et al. (1995 and in preparation). 
The elliptical data (Lucey & Carter 1988; Faber et al. 1989; 
Dressier, Faber & Burstein 1991) are the same data previ- 
ously compiled by Burstein (1990) for the Mark II catalogue. 
The spiral data contain the nearby Mark II data of Aaronson 
et al. (1982) as revised by Tormen & Burstein (1995), and 
newer extended data from Han (1991), Han & Mould (1992), 
Mould et al. (1991, 1993), Willick (1991), Courteau (1992), 
and Mathewson, Ford & Buchorn (1992). The sample con- 
sists of about 3000 galaxies grouped into ~1150 objects. 

Inhomogeneous Malmquist (IM) bias affects the deter- 
mination of Q (or P) from TF data when the analysis, 
as in potent, combines neighbouring galaxies in inferred- 
distance space (e.g. Willick 1991, 1994; Hudson 1994a; Dekel 
1994). The distance errors, combined with galaxy-density 
variations along the line of sight, systematically enhance 
the inferred velocity gradients and thus the inferred den- 
sity fluctuations and SI or /3. The IM bias is first reduced 
and then corrected by a new procedure that has been tested 
using mock peculiar velocity catalogues drawn from N-body 
simulations. The mock catalogues are constructed in a way 
which mimics the non-uniform sampling and the noise of 
the real data. The IM bias correction procedure is discussed 
in detail by Dekel et al. (in preparation, hereafter D95c) 
and is only summarized briefly here. First, the galaxies are 



grouped heavily in redshift space, while correcting for 'se- 
lection' bias (Willick 1991, 1994). This grouping reduces the 
distance error of each group of N members to A/\/iV, where 
A ~ 0.15 — 0.21 and 0.21 for spirals and ellipticals respec- 
tively. The IM bias scales like the square of the distance er- 
ror, so this grouping significantly reduces the IM bias. Then, 
the noisy inferred distance of each object, d, is replaced by 
the expectation value of the true distance r given d (e.g. 
Willick 1991, equation 5.70): 



E(r\d) 



f o °° r 2 n(r) exp {- [ln( ^f ]2 ^j dr 



(2) 



For single, ungrouped galaxies, n(r) should be the number 
density of galaxies in the underlying distribution from which 
galaxies were selected for the sample (assuming they were 
selected by quantities that do not explicitly depend on r). In 
cases where the sample is redshift-limited, n(r) is truncated 
at the appropriate distance. In a grouped sample, the den- 
sity run n(r) is multiplied by appropriate grouping/selection 
correction factors for the grouped and ungrouped objects re- 
spectively, (see D95c for details). This procedure is found to 
reduce the IM bias in the mock catalogues to the level of a 
few per cent. As an approximation to n(r) we use for spi- 
rals the density of IRAS 1.2-Jy galaxies (Fisher 1992), and 
for ellipticals the density of early-type galaxies derived by 
Hudson (1993a, see also Section ^| below), both smoothed 
with a Gaussian filter of radius 500 km s . The final IM 
correction typically amounts to less than 10 per cent in the 
1200 km s _1 smoothed density field described next. 

The potent procedure takes these discrete and noisy 
radial velocity data, u, at positions Xi, and first smoothes 
them with a Gaussian window of radius 1200 km s _1 (DBF; 
Dekel et al., in preparation, hereafter D95a, D95b). It then 
applies the gravitational ansatz of potential flow to re- 
cover the missing transverse components of the velocity field 
(Bertschinger & Dekel 1989). Finally, it applies a quasi- 
linear approximation with an assumed value of SI to re- 
construct the associated field of mass-density fluctuations. 
The quasi-linear generalization of the linear approximation 
to gravitational instability is the solution of the continu- 
ity equation under the Zel'dovich approximation in Eulerian 
space (Nusser et al. 1991), 



5 c (x) = ||I- f^dv/dxW - 1 



(3) 



where the bars denote the Jacobian determinant, and I is 
the unit matrix. 

The non-trivial step in potent is the smoothing of the 
data into a radial velocity field u(x). This smoothing proce- 
dure is described and tested in detail elsewhere (DBF; D95b; 
Kolatt et al., in preparation). The procedure is only briefly 
summarized here with emphasis on the new features of the 
method. The aim is to reproduce the u(x) that would have 
been obtained had the true three-dimensional velocity field 
v(x) been sampled densely and uniformly and smoothed 
with a spherical Gaussian window of radius i? s . With the 
discrete data available, u(x c ) is taken to be the value at 
x = x c of an appropriate parametric local velocity model 
v(ak, x~ x c ) obtained by minimizing the weighted sum of 
residuals 
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Wi [Ui - Xi ■ v(ce k ,Xi — x c 



(4) 



in terms of the parameters ah within an appropriate local 
window Wi — W(xi,x c ). The local velocity model and the 
weighting needed in order to minimize random errors and 
systematic biases are chosen as follows. 

Even for the case of dense, uniform sampling the re- 
covered u(x) field will in general suffer from 'window bias'. 
Unless Rs is much smaller than the distance from the origin 
to the centre of the window, r c , the radial peculiar velocity 
data, u iy cannot be averaged as scalars because the direc- 
tions Xi differ from x c , so tt(ir c ) requires a fit of a local 3D 
velocity model. The original potent method (DBF) used 
the simplest local model, v(x) = B, where B is a uniform 
velocity, for which the solution can be expressed explicitly in 
terms of a tensor window function. The tensorial correction 
to the spherical window has conical symmetry, weighting 
more heavily objects of large Xi-x c . This window deforma- 
tion introduces a bias which is particularly severe nearby, 
where R s /r c is not small, and in places where the velocity 
has a large divergence or convergence transverse to the line 
of sight. For example, a converging flow in a plane transverse 
to the line of sight would be biased to show an artificial com- 
ponent towards the origin. For the current peculiar velocity 
data, the resulting bias in the smoothed radial peculiar ve- 
locity field reaches a maximum value of ~ 300 km s _1 at 
the position of the Great Attractor. A way to reduce this 
bias is by generalizing v(x) into a linear velocity model, 
v(x) = B+~L-(x— x c ), with L a symmetric tensor, which en- 
sures local irrotationality. The linear terms tend to 'absorb' 
most of the bias, leaving v(x c ) = B less biased. Unfortu- 
nately, a high-order model tends to pick undesired small- 
scale noise, especially where the data are sparse and noisy. 
The optimal compromise is found to be a first-order model 
fit out to r = 40/i _1 Mpc, smoothly changing to a zeroth- 
order fit beyond 6O/1" 1 Mpc (D95b). 

Unfortunately, the real peculiar velocity data are both 
non-uniformly sampled and noisy. If the true velocity field is 
varying within the effective window, then non-uniform sam- 
pling will introduce a sampling-gradient (SG) bias. This bias 
occurs because the smoothing is weighted by the distribu- 
tion of sampled galaxies whereas the aim is equal-volume 
weighting. One should weight each object by the local vol- 
ume it 'occupies', Vi. We adopt Vi oc R\ where Ra is the 
distance to the fourth neighbouring object. This weighting 
procedure (termed P0) is found from simulations to reduce 
the SG bias in the Mark III sample to negligible levels out 
to a typical depth of 60fo -1 Mpc (away from the Galactic 
Zone of Avoidance). The Ri(x) field will serve later as a 
flag for poorly sampled regions, to be excluded from any 
quantitative analysis. 

The ideal weighting for reducing the effect of Gaus- 
sian noise in the peculiar velocity measurements, at, requires 
weights Wi(xa~ 2 , but this spoils the volume weighting and 
biases u towards its values at smaller Vi and at nearby clus- 
ters where the errors are small. A successful compromise 
(termed PI) is to weight by both, i.e. 



W(xi,x c ) oc Vicr^ 2 exp[— (xi — x c ) 2 /2R 2 



(•») 



magnitude of the systematic uncertainty still remaining in 
the potent recovery. 

The errors in the recovered fields, due to scatter in the 
distance indicator and measurement errors, are assessed by 
Monte Carlo simulations, where the input distances are per- 
turbed at random using a Gaussian of standard deviation 
Oi [in fact (of + a 2 ) 1//2 , where 01 ~ 200 km s _1 is the local 
dispersion of velocities in the 'field'] before being fed into 
potent. The error in 8p at a grid point is estimated by the 
standard deviation of the recovered <5p over the Monte Carlo 
simulations, as- In the well-sampled regions, which extend 
in Mark III out to 40-60h~ 1 Mpc, the errors are a s ^0.1-0.3, 
but they may blow up in certain regions at large distances. 
To exclude noisy regions, any quantitative analysis should 
be limi ted to points where as is within certain limits (see 
Section 



4.1 



below) . 

There exists a fundamental freedom in determining the 
zero-point TF parameter which fixes the distance scale in 
units of km s _1 . A change of this zero-point multiplies all 
distances by a factor (1+e) while the redshifts are fixed. This 
is equivalent to adding a monopole Hubble-like component 
— er to the peculiar velocities v, and an offset 3e/ _1 (f2) to 
S (~ — /~ 1 V ■ v). The zero-point Hubble flow has tenta- 
tively been determined in the Mark III data from the 38 
spiral clusters, which are distributed across the the sample 
volume. After processing through potent, the zero-point is 
fine-tuned by minimizing the volume-weighted variance of 
the recovered peculiar velocity field (equivalent in the lin- 
ear regime to enforcing (S) — 0) in a 'fair' volume, taken 
here to be a sphere of radius 6000 km s _ . In fact, the zero- 
point ad opte d has no effect on the determination of /3 below 
(Section O). 



The difference between the recovered fields using the P0 
and PI weighting schemes can serve as an indicator for the 



The potent density field (D95a) is shown in the left- 
hand panels of Figs |l| and |^. Fig. jj] shows contour plots of the 
smoothed density field in slices parallel to the Supergalactic 
Plane. Fig. |^ shows the density field in the Supergalactic 
Plane as a surface plot where the height of the surface above 
the plane of the plot is proportional to density. 



3 THE OPTICAL DENSITY FIELD 

We now summarize the optical galaxy sample, complete- 
ness corrections and the method used to reconstruct the 
density field of optical galaxies (for further details see Hud- 
son 1993a). The basic catalogue of 12 439 galaxies covering 
67 per cent of the sky is a merger of the diameter-limited 
UGC (Uppsala General Catalogue of Galaxies: Nilson 1973) 
and ESO (The ESO-Uppsala Survey of the ESO(B) Atlas: 
Lauberts 1982) catalogues. The UGC catalogue is corrected 
for incompleteness for diameters between 1.2 and 1.6 arcmin, 
and the ESO catalogue is complete at diameters larger than 
1.4 arcmin (Hudson & Lynden-Bell 1991). The 6747 redshifts 
in the sample were obtained from the major redshift sur- 
veys, i.e., CfAl (Davis & Huchra 1982; Huchra et al. 1983), 
SSRS (da Costa et al. 1988, 1991), SPS (Dressier 1991), and 
Perseus-Pisces (Giovanelli & Haynes 1985; Giovanelli et al. 
1986), and from the redshift compilations of Huchra (1990) 
and Fairall & Jones (1991). The variable completeness of the 
redshift sample is mapped as a function of angular diameter 
and position on the sky. These completeness functions are 
used in conjunction with the UGC and ESO diameter func- 
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Figure 1. The density fields in slices parallel to the Supergalactic Plane. The left-hand panels show the POTENT mass density field 
and the right-hand panels show the density field of optical galaxies, both smoothed with a Gaussian filter of radius 1200 km s —1 . The 
contour levels are in steps of A<5 = 0.2. The medium-thickness contours indicate the mean density, thin solid contours show overdense 
regions and dashed contours show underdense regions. The heavy contours indicate the boundary of the 'Large' comparison volume (see 
Section for details'). The dotted circle has a radius of 8000 km s 1 . The toD row shows the slice at SuDcrealactic Z = 2000 km s 1 . 
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Figure 2. Surface plots of the density fields in the Supergalactic Plane. The left-hand panel shows the POTENT mass density field and 
the right-hand panel shows the density field of optical galaxies, both smoothed with a Gaussian filter of radius 1200 km s~ x . The density 
contrast is proportional to the height of the surface above (or below) the plane of the plot. 



tions (Hudson & Lynden-Bell 1991) to assign a weight to 
each galaxy in the redshift sample to compensate for those 
galaxies in the underlying density field that were not se- 
lected. 

We then iteratively transform from redshift space to 
distance space by self-consistently calculating linear peculiar 
velocities from the optical density field (cf. Yahil et al. 1991). 
The value of fio is fixed at 0.5 for these iterations, but the 
resulting optical density field does not depend sensitively on 
this choice. Indeed, had we used f3o = 0.75 in the iteration, 
our final result for /3o would have increased by only 0.03. 

There are no data in the Zone of Avoidance |6| < 12° 
and the 'missing' equatorial strip between the UGC and 
ESO catalogues (-17? 5 < Dec. < -2? 5). We adopt a simple 
scheme for filling in this 'masked' area: the 'cloned mask' 
model (cf. Lynden-Bell, Lahav & Burstein 1989) attempts 
a crude interpolation in the galactic and equatorial miss- 
ing strips by copying the set of weighted galaxies in adja- 
cent strips and shifting these 'clones' in galactic latitude and 
declination respectively. 

In the final step, the weighted set of points is smoothed 
to yield an optical density field. We use the same smooth- 
ing filter as used for the potent density field, namely a 
Gaussian filter with a 1200 km s _1 radius, and we sample 
the density fields on the same 500 km s _1 grid. In the final 
comparison, we excl ude all grid points within the masked 
area (see Section 4.1 below). However, note that, due to the 
smoothing, the density assumed within the masked region 
will affect the density at adjacent grid points. Therefore, the 
'cloning' procedure is preferred over less realistic schemes 
such as the 'average mask model' in which the masked re- 
gion is set to the average density. 

The mean density of optical galaxies is needed in order 
to calculate the density contrast, So — (po — po)/pcT- I n this 
paper, we estimate the mean optical density by averaging 
over optical galaxies in the unmasked volume within 6000 
km s _1 . This volume is close to t he t ypical depth of the final 
comparison volume (see Section 4.1 below). 

The optical density field is shown in the right-hand pan- 
els of Figs [l] and |^. 



4 DETERMINATION OF f3 

In the determination of the parameter fio, our philosophy 
is to carry out as simple an analysis as possible. We use 
simplifying approximations where we believe, based on our 
experience with mock data sets, that any biases introduced 
by these approximations are small compared to the unavoid- 
able uncertainties due to the random errors, the limited vol- 
ume sampled (cosmic scatter), the unknown complexity of 
the true biasing relation, and the remaining uncorrected ef- 
fects of non-linear gravity. We deviate from a straightfor- 
ward analysis only when faced with a significant bias. The 
approximations made in this analysis are tested using mock 
catalogues elsewhere (Dekel et al., in preparation, hereafter 
D95d). 

we define the volume used for the com- 
In Section 4.2, we discuss the method of fit, based 
on a simple \ 2 statistic, and the effective number of inde pen- 



In Section |4.1 
parison 



dent volumes within the comparison volume. In Section 4.3 



we discuss the issue of SG bias and how it biases the deter- 



mination of Po, and in Section 4.4 w e di scuss the scatter in 
the density-density plots. In Section 4.5, we present results 
for Po assuming fl = 1, and in Section 4.6 we discuss the 
results for SI = 0.3. 



4.1 Definition of the comparison volume 

The density-density comparison in this paper is local (unlike 
the velocity-velocity comparison) and the potent data are 
weighted inversely by their local errors. Consequently, we are 
free to choose the volume considered in the density-density 
comparison. In selecting the volume, we aim to minimize the 
SG bias and the random errors, and to maximize the number 
of effectively independent volumes used in the comparison. 
In practice, there is an obvious trade-off between these two 
goals. We therefore limit the volumes by R4, in order to 
minimize the SG bias, and by as, in order to eliminate re- 
gions with very large random errors. Furthermore, we use 
only the grid points that are within 7000 km s _1 and out- 
side the masked region of the optical density field. In order 
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Table 1. 


Comparison 


volumes. 






Sample 


R 4 






N eS 




km s" 1 








PI S 


< 1200 


< 0.25 


2197 


22.1 


PI L 


< 1500 


< 0.30 


3334 


31.2 


PO S 


< 1200 


< 0.30 


1712 


18.0 


PO L 


< 1500 


< 0.40 


2934 


28.9 



S = iV, 



to investigate how our results depend on specific choices for 
these parameters, we consider two comparison volumes for 
each of the potent methods (P0 and PI). We refer to these 
as the S ('small') and L ('large') volumes. Table |l| lists the 
7?4 and as limits, the number of grid points (from a cubic 
grid of 500 km s _1 spacing) used in the comparison, and the 
effective number of independent points in the fit (discussed 
below) for the four volumes considered. 

The limits of the P1L volume are indicated in Fig. |l| 
The radial indentations in the comparison volume arise from 
the optical mask (i.e. the Zone of Avoidance and the 'miss- 
ing' equatorial strip) . The irregular outer limit is determined 
by the inhomogeneity in the potent data. The PIS volume 
covers approximately 67 per cent of the sky and has a total 
volume of (6500 km s -1 ) 3 . Its outer boundary (outside of 
the Zone of Avoidance) lies between 3000 and 7000 km s , 
depending on direction. 



4.2 Method of fit 

A visual comparison between the left- and right-hand panels 
of Fig. |l| and Fig. ^| shows a reasonable resemblance between 
the large-scale features of the potent and optical density 
fields. Displacements of density peaks are seen at a similar 
level in the mock catalogues as a result of the sparse, non- 
uniform and noisy sampling. Assuming GI and no systematic 
errors, the mass density - optical density diagram is a direct 
representation of the biasing relation (for the value of Q used 
in the potent reconstruction). Assuming a linear biasing 
relation with no scatter, we expect the densities to be related 
by 



8p = X5o + c 



(6) 



where the slope, A, is just equal to &Q 1 . The zero-point, 
c, allows for a relative offset in the mean density of the 
two density fields. This parameter is necessary because the 
volumes within which we have separately normalized the 
potent and optical density fields are slightly different from 
each other. 

The random errors in the inferred peculiar velocities in- 
troduce an error in the potent density, as, so of course we 
expect some scatter in the potent - optical density dia- 
gram. We have estimated the random shot-noise errors in 
the optical density field using the technique of bootstrap re- 
sampling. We find these errors to be typically ~ 0.07, and 
always less than half the random errors of the potent den- 
sity field at a given grid point. In the following analysis we 
assume that all errors are in the potent density field, and 
thus regress 8p on <5o in the determination of A. 

If the grid points were independent of each other and if 
potent errors were Gaussian, we could construct a reduced 
\ 2 statistic, 



grid 



E 



(Sp - XSo - cf 



(7) 



where the sum is over all iV g rid grid points in the comparison 
volume and ap is the error in the potent density. 

In practice, we sample the density fields on a grid with a 
much finer spacing (500 km s _1 ) than the smoothing radius 
(1200 km s _1 ), so the grid points are not independent. The 
effective number of independent volumes in the \ 2 fit can 
be estimated by 



N-, 



N. 



si id 



EE 



exp( 



(8) 



(DBY93). The values of N e g for each comparison volume 
are given in Table [j]. We find that with our grid spac- 
ing and smoothing we have oversampled the volume by 
a factor F = N gI i<i/N e ff ~ 100. Thus our data points 
would be independent if sampled on a grid with spacing 
F 1/3 x 500 km s" 1 = 2320 km s"\ which is about twice 
the smoothing length. Due to the oversampling, the appro- 
priate statistic to use is a corrected XeS = N e ffS, which is 
distributed like \ 2 with Nd.o.t = N c « — 2 degrees of free- 
dom. This statistic is used to assess the errors on A and c, 
as a simple approximation to the elaborate and more accu- 
rate likelihood analysis using Monte Carlo simulations (e.g. 
DBY93). 



4.3 Correction of biases in potent 

Fig. ^ shows the comparison between Sp (PIS) and So- If 
we set ap = as (the random error in the potent density) 
in equation (Q), we obtain the result A = 0.66 ± 0.08, c = 
—0.05 ± 0.02. The quoted errors are the formal errors from 
the x 2 fit- This fit is shown as the solid line in Fig. || 

However, it would be premature to adopt these results 
as our best estimate of f3o because this estimate is known to 
be biased. DBY93 tested the potent method by comparing 
the smoothed IRAS density field to the density field recov- 
ered by the potent method from the unsmoothed IRAS- 
predicted peculiar velocities with Q — 1 and bi = 1. The 
IRAS velocity and density fields served there as a model for 
a general gravitating system. DBY93 found that, in the limit 
of dense, uniform sampling, they recovered a slope close to 
the true slope of unity with very small scatter (see e.g. their 
fig. 4c). However, when they sampled the IRAS peculiar ve- 
locity field at only the positions of the galaxies in the Mark 
II peculiar velocity data sample, they found an increased 
scatter and a biased slope of 0.66. Thus, the systematic er- 
rors in the potent procedure are correlated with the density 
field and introduce a bias in the slope, tending to reduce it. 

Following DBY93, we use as a test case the optical den- 
sity field itself from which we calculate predicted peculiar ve- 
locities via gravity at the estimated distance of each galaxy 
in the peculiar velocity data set. Specifically, we begin by 
smoothing the optical density field with a Gaussian filter of 
500 km s" 1 radius. Peculiar velocities are calculated using 
the prescription of Nusser et al. (1991), i.e. 



>(r) = 



MO 



4tt&o J 1 + O.lSfc^Mr') \r' -r 



3d r , 



(9) 
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Figure 3. The density-density scatterplot with the optical density (O) on the horizontal axis and the POTENT mass density (P) on 
the vertical axis. The densities are sampled on a 1000 km s _1 grid within the 'PI Small' volume, thus the plotted points oversample 
the number of independent volumes by a factor ~ 12. The areas of the plotted points have been scaled with the inverse squares of their 
errors. The solid line is the result of the weighted linear regression of <5p on 8q; a dotted line of slope unity is also shown for reference. 



at the observed positions of the galaxies in the Mark III 
data set. For the = 1 case, we adopt a biasing param- 
eter of bo = 1.3, which ends up consistent with our final 
result. These predicted peculiar velocities are then used as 
the input to the potent procedure with 1090 km s _1 Gaus- 
sian smoothing which, when added in quadrature with the 
500km s" 1 smoothing radius of the input field, ensures a to- 
tal Gaussian smoothing of 1200 km s . The resulting mass 
density fluctuation field is then multiplied by bo to obtain 
the corresponding field of optical galaxies. The result of this 
procedure, which we label P(O), is a galaxy density field 
which is affected by SG bias in the same way as the potent 
mass density field recovered from the Mark III data. Fig. ^ 
compares the density fields of potent [P], optical through 
potent [P(O)], and raw optical [O]. 

In Fig. H we show the P(O) - O scatterplot for this 
test case. If the same weights at each grid point as used 
in the P - O comparison above are adopted, the slope of 
best fit is 0.88. This bias in the slope is smaller than that 
found by DBY93 due to the increase in sampling density of 
the peculiar velocity data, the improvements in the potent 



procedure, and the more accurate adjustment of the potent 
smoothing scale to a total of 1200 km s -1 . 

4.4 Scatter in the potent — optical density 
diagram 

Before turning to the results of the fits, we must first dis- 
cuss the scatter in the density-density diagram, in order to 
verify that the assumed model of GI plus linear biasing is 
indeed consistent with the data. We find that, by setting 
op equal to the error in the Monte Carlo simulations ag, 
we obtain \eS ii 2iV c ff , which is apparently too large for an 
acceptable fit. This could be a result of underestimating the 
errors. Indeed, it is likely that there are other sources of er- 
ror, predominantly in the POTENT density field, which are 
not accounted for by the Monte Carlo perturbations of the 
distances. We therefore allow for additional errors via an ad- 
ditional parameter, <7 a dd- The additional errors may or may 
not be correlated with as, so we assume a weak correlation 
by adding linearly: ap = er ad d + o*. We then adjust a ad d so 
that Xeff ~ ^cff — 2. The resulting values are <r ad d = 0.06 
for the PI comparisons, and <r a dd = 0.04 for the P0 compar- 




Figure 4. Maps of the density field in the Supergalactic Plane. The contour levels are as in Fig. 1. The upper left-hand panel shows 
the POTENT mass density field (P). The upper right-hand panel shows P(O), the POTENT recovery using predicted peculiar velocities 
from the optical density field at the positions of the Mark III data. The lower right-hand panel shows the optical density field (O). The 
difference between P(O) and O shows the effect of the biases discussed in the text. 



isons. In Section 4.5, we show that the value of the slope, A, 



is not significantly affected by our choice of cr a dd- It is none 
the less necessary to include this additional error in order 
to obtain errors on the slope that are consistent with the 
observed scatter in the density-density plot. 

Are these values of a a dd reasonable compared to the 
possible sources of scatter that are not included in as, and 



which are not eliminated by the comparison of P with P(O)? 
Additional sources of error might include the following. 

(i) The Monte Carlo estimates, as, may underestimate 
the total errors in the potent density field. We have made 
a preliminary test of these errors using an ungrouped version 
of the mock catalogue described above. We have compared 
densities obtained from applying the potent procedure to 
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Figure 5. The density— density scatterplot with optical density (O) plotted on the horizontal axis and P(O), the POTENT recovery 
using predicted peculiar velocities from the optical density field at the positions of the Mark III data, plotted on the vertical axis. The 
comparison volume, sampling and symbols are as in Fig. 3. The solid line is the weighted linear regression of <5p(o) 011 ^Oi a dotted line 
of slope unity is also shown for reference. The deviation from a slope of unity is due to the biases discussed in the text. 



the unperturbed velocities of galaxies at their true distances 
and to the perturbed velocities at the perturbed positions, 
after correcting for IM bias. We find that in addition to the 
Monte Carlo errors, as, we require <r a dd = 0.03 in order to 
obtain an acceptable value of \ 2 - 

(ii) There may be errors in the potent density field due 
to imperfections in the inhomogeneous Malmquist bias cor- 
rections. An upper limit to these errors is estimated by the 
rms difference between the potent density fields corrected 
and uncorrected for IM, which is ~ 0.11. 

(iii) There may be systematic errors in the optical den- 
sity field due to interpolation ('cloning') into the unsampled 
regions. An upper limit to these errors is estimated by com- 
paring the cloned mask density field to the average mask 
density field. Within the PIS comparison volume, the rms 
scatter between the cloned and average mask density fields 
is ~ 0.06. Spurious large-scale gradients may also exist at 
some level due to possible errors in matching the UGC and 
ESO data in the northern and southern hemispheres (see 
Hudson 1993b, 1994b). We note, however, that the results 
for Po from these two hemispheres sepa rate ly are in perfect 
agreement with each other (see Section 4.5 below). 



Table 2. Results of P - P(0) fits. 



Sample 

PIS 

P1L 

P0S 

POL 



A 

0.74±0.13 

0.78±0.13 
0.65±0.15 
0.68±0.14 



-0.01±0.04 

-0.02±0.04 
0.01±0.04 
0.01±0.04 



(iv) The shot-noise error in the optical density, although 
smaller than the random errors in potent, is not negligible. 
Within 6000 km s , the typical shot-noise error is ~ 0.07. 

(v) Spatial variations in the efficiency of galaxy formation 
and deviations from a linear biasing relation would generate 
true, physical, scatter about an assumed linear biasing rela- 
tion between galaxies and mass. A scatter of ~ 0.15 in the 
optical galaxy density added in quadrature to the potent 
error would account for the additional scatter in Fig. |j| 

Since these possible sources of scatter are of the same 
order as our adopted value of <T a dd, we may conclude that 
there is no evidence for inconsistency between the data and 
the assumptions of GI and linear biasing, allowing us to 
proceed to the parameter fit. 
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4.5 Results 

The comparison between the potent density and the 
POTENT-recovered optical density P(O) removes the biases 
in the potent procedure and thus yields an unbiased esti- 
mate of /3o- As discussed above, we assume that the errors 
in the potent density field dominate. We therefore obtain 
the slope from a regression of Sp on <5p(o) with the potent 
error op = os + <7 a dd. The results of the fit for the PO and PI 
(L and S) cases are summarized in Table ^. Fig. ^| shows the 
P - P(O) scatterplot for the PIS case. In all cases, the zero- 
point c is small and compatible with zero. We have made 
a number of tests to check the robustness of the value of A 
derived from the regression of 5p on <5p(o)- The results of 
these tests are summarized below. 

(i) The value of A is insensitive to the details of the as- 
sumed POTENT errors, provided that the optical errors are 
negligible. For the PIS case, if we assume equal errors for 
all grid points, our resultant A increases by only 3 per cent. 
Similarly, if we set cr a dd = 0, then A decreases by 5 per cent. 
If we fit by minimizing absolute deviation, which is less sen- 
sitive to outliers, then A decreases by 3 per cent. 

(ii) The value of A is insensitive to the exact volume 
used, as indicated by the good agreement between the S 
and L volumes. A further check can be made by dividing 
the PIS volume into the northern celestial hemisphere (con- 
taining the Perseus-Pisces supercluster) and the southern 
hemisphere (containing the Great Attractor region) . We find 
A = 0.73 ± 0.26 and 0.75 ± 0.15 for northern and southern 
hemispheres, respectively. In fact, this almost perfect agree- 
ment between these independent volumes may indicate that 
we have overestimated our random errors. 

(iii) The results are insensitive to the details of the 
Malmquist corrections. As a test of these corrections, we 
have compared the optical density field to a potent density 
field in which the input data were corrected for Malmquist 
bias assuming a homogeneous galaxy density field (here- 
after HM-corrected). Although in a few specific regions (e.g. 
Perseus-Pisces) the difference between the IM- and HM- 
corrected potent density fields is as large as ~ 30 per cent 
in Sp, when averaged over all grid points in the compari- 
son volume the net effect is small. We find that with the 
HM corrections A is only 2 per cent higher than when IM 
corrections are used. 

(iv) The values for A are ~ 12 per cent lower for the P0 
case than for the PI case. This is a crude measure of the 
typical systematic error associated with the potent proce- 
dure. 

(v) We have compared the PIS case to a version of the 
optical density field in which the masked region was set to 
the average density before smoothing (the 'average mask 
model' from Hudson 1993a). The value of A is ~ 15 per cent 
higher than the result with the 'cloned mask' density field 
used here. Although the average mask model is unrealistic, 
it sets an upper limit on the systematic error arising from 
the treatment of the masked region. 

Our estimate of j3o, adopting the PIS results of the 5p 
on <5 P (o) regression, is fio — 0.74±0.13, where the quoted la- 
errors are the formal errors on the slope which are consistent 
with the scatter about the linear fit. 



Note that if we assume that the 'extra' scatter in the 
density-density diagram lies in the optical density field, 
whether due to errors in the optical density field or due 
to a real physical effect, such as spatial variations in the ef- 
ficiency of galaxy formation, then we must do a more elab- 
orate regression which allows for errors in both the potent 
and optical density fields. As an example, if we set er a dd = 
and thus op = os, but allow for a 0.15 error in the optical 
density, the regression yields a slope A = 1.00 ± 0.17 and an 
acceptable value of XeS- ^ n Section |^ we discuss the possi- 
bility that there may be significant scatter in the efficiency 
of galaxy formation. 

4.6 The dependence on Q 

The potent reconstruction includes quasi-linear effects so, 
in principle, it should be possible to obtain some information 
on Q, independent of bo- To test this idea, we have generated 
a potent reconstruction with SI set to the value 0.3. This 
is compared to a P(O) density field in which the predicted 
peculiar velocities are determined via equation (^) with SI = 
0.3, and bo = 0.7 in the quasi-linear correction within the 
integral. The slope of the resulting fit is thus an estimate of 
b Q \ 

We find a slope A = 6 D = 1.63 ± 0.33, corresponding 
to bo = 0.61 ± 0.12 or (3 = 0.79 ± 0.16. Thus, we recover 
essentially the same value of j3o as in the SI = 1 case because 
we are only in the mildly non-linear regime. The scatter 
about the straight-line fit in the SI = 0.3 case is marginally 
smaller than in the SI = 1 case, but the improvement is not 
statistically significant. We conclude that, given the random 
and systematic errors, including the possible deviations from 
linear biasing, we cannot remove the degeneracy between SI 
and bo- In order to rule out a low-density universe by these 
results, one has to appeal to external arguments such as the 
fact that the required degree of anti-biasing is physically 
implausible. 



5 DISCUSSION AND SUMMARY 

Our simple analysis does not allow a formal evaluation of 
goodness of fit between the data and the assumed model of 
GI and linear biasing, but we conclude, based on the ex- 
pected random and systematic errors and the possible real 
scatter in the biasing relation, that the data and model are 
not in conflict. Note, however, that a non-gravitational ori- 
gin to the peculiar velocities, which could invalidate the cur- 
rent analysis, cannot be ruled out. An apparent agreement 
between the potent density field and the galaxy density 
field can also occur in non-gravitational scenarios provided 
that the galaxy density field has developed from homoge- 
neous initial conditions in a way that satisfies the continuity 
equation, and that the present-day velocity field is irrota- 
tional and proportional to the time-averaged velocity field 
(Babul et al. 1994). 

Adopting the model of GI and linear biasing, our main 
goal has been to determine (5o on a scale of 1200 km s" 1 . We 
find, from a regression of potent density on optical density, 
the result j3o = 0.74 ± 0.13 (formal lo error). This result is 
robust to the details of the volume used, to the treatment of 
Malmquist bias, and to the actual value of SI. The regression 
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Figure 6. The density— density scatterplot with P(O) plotted on the horizontal axis and the POTENT density (P) plotted on the vertical 
axis. The comparison volume, sampling and symbols are as in Fig. 3. The solid line is the straight line of best fit; a dotted line of slope 
unity is also shown for reference. The slope A is our estimate of b^ 1 for Q = 1. 



of potent density on optical density is insensitive to the 
relative weighting of grid points in the fit. The sensitivity of 
(5o to extreme changes in the potent weighting scheme is 
at the 12 per cent level. Its sensitivity to extreme changes 
in the way in which the optical density field is interpolated 
across of the Zone of Avoidance is at a similar level. 

Non-linear effects on this scale are small compared to 
the uncertainties in the data, the analysis, and the biasing 
scheme, so we are unable to remove the degeneracy between 
ft and &o. Our degenerate result corresponds to bo — 1.35 ± 
0.23 if ft = 1, and to ft = 0.61 ± 0.18 if bo = 1. Low values 
of the density parameter (ft ^ 0.3) would require optical 
galaxies to be significantly anti-biased with respect to mass. 

One should be cautioned not to interpret the estimated 
value of /3o too literally. On top of the formal error of 0.13 
and the systematic uncertainties of similar order, our esti- 
mate of Po is subject to cosmic scatter in the mean density; 
the comparison volume of radius ^ 6000 km s _1 may still 
be small for a 'fair' sample. This cosmic scatter arises be- 
cause, in our analysis, both the peculiar velocity data input 
to potent and the mean density of optical galaxies are nor- 
malized within a volume of radius ~ 6000 km s . Note that 
while the free parameter c corrects for small differences in 



the normalizing volumes of the potent and optical den- 
sity fields, and consequent differences in the mean density 
of these fields with respect to each other, it does not nec- 
essarily correct to the universal density. Thus our result is 
a local value of Q measured within the comparison volume. 
The rms scatter of the mean density in such a volume is 
expected to be at the level of ~ 10 — 15 per cent according 
to the fluctuation power spectra in conventional cosmologi- 
cal scenarios. Thus the local value of ft may differ from the 
universal value at a similar level. 

We have assumed throughout that the additional scat- 
ter in the density-density scatterplots (over and above that 
due to the distance errors which scatter the potent den- 
sities) is due to other small unquantified e rror s in the PO- 
TENT density fields. As discussed in Section 4.4, the amount 



of additional scatter is comparable to our estimates of these 
systematic errors. However, it is also possible that spatial 
variations in the biasing relation between galaxies and mass 
may be large on scales of 12h~^ Mpc. For example, Cen 
& Ostriker (1993) investigated hydrodynamic simulations of 
the Cold Dark Matter cosmology and found on scales of 
10ft" 1 Mpc a non-linear biasing relation which corresponds 
roughly to bo — 1-2 — 1.4 in the range —0.5 < 8 < 1.0, in 
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good agreement with our results. The rms scatter in their 
biasing relation is about 0.23 in So at So = 0, compara- 
ble to the scatter of ~ 0.15 which we would need to add in 
quadrature to the optical galaxy density to obtain good fits. 
If there are significant spatial variations in the efficiency of 
galaxy formation, then this adds an astrophysical source of 
scatter in the 5 - So plot which must be accounted for in 
the analysis. If this extra scatter is assigned to <5o then the 
slope of the regression changes to /3o ~ 1. 

The current estimates of 57 and b from large-scale struc- 
ture are reviewed by Dekel (1994, section 8 and his table 1). 
In particular, the result for f3o found here is lower than 
the value Pi = 1.28+o'gg (95 per cent confidence) found by 
DBY93, who compared an earlier realization of potent with 
the 1.9-Jy IRAS density field. A comparison between the 
1.2-Jy IRAS density field and the new version of potent 
is in progress (D95d). If we allow for the fact that IRAS 
galaxies are less clustered (biased) than optical galaxies by 
a relative factor of ~ 1.3 on the scales studied in this paper 
(Hudson 1993a; Peacock & Dodds 1994), then the different 
results for P lead to compatible values of Q. 

Hudson (1994b) compared predicted peculiar velocities 
in linear theory from the optical density field smoothed 
with a 500 km s _1 top-hat filter to observed peculiar ve- 
locities from the Mark II catalogue, and obtained the result 
Po = 0.50 ± 0.06 (where the quoted errors are formal la- 
random errors only). It seems at this early stage, while the 
various analyses are still in progress and the systematic ef- 
fects are not fully understood, that velocity-velocity com- 
parisons tend in general to yield somewhat lower estimates 
for (3 than the estimates based on density-density compar- 
isons (e.g. Nusser & Davis, private communication; Willick 
et al., in preparation). However, the results are not at all in- 
compatible when one allows for the systematic uncertainties 
in both determinations. A possible physical explanation for 
this difference in /3 is a variation in the biasing relation with 
scale or density. 

The lower bounds on Q from cosmic flows independent 
of biasing (e.g. Nusser & Dekel 1993; Dekel & Rees 1994; 
see the review by Dekel 1994) are currently at the level of 
Q > 0.3 with 99 per cent confidence, consistent with the the- 
oretically favoured value of 57 = 1. Other input of relevance 
include (i) constraints from the microwave background fluc- 
tuations, (ii) rms number fluctuations of optical and IRAS 
galaxies in top-hat spheres of radius 8/i -1 Mpc (as,o ~ 0.95, 
&8,i ~ 0.65), and (iii) theoretical expectations for biasing 
based on first-generation simulations with gas dynamics. A 
choice of parameters that is consistent with the result of 
this paper and all the constraints above would have mass 
fluctuations at the level of as « 0.75, and linear biasing fac- 
tors bo « 1.3 and bi w 0.9 for optical and IRAS galaxies 
respectively. 
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